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Abstract. We investigate the influence of aperiodic modulations of the exchange interactions between 
nearest-neighbour rows on the phase transition of the two-dimensional eight-state Potts model. The sys- 
tems are studied numerically through intensive Monte Carlo simulations using the Swendsen-Wang cluster 
algorithm for different aperiodic sequences. The transition point is located through duality relations, and 
the critical behaviour is investigated using FSS techniques at criticality. While the pure system exhibits 
a first-order transition, we show that the deterministic fluctuations resulting from the aperiodic coupling 
distribution are liable to modify drastically the physical properties in the neighbourhood of the transition 
point. For strong enough fluctuations of the sequence under consideration, a second-order phase transition 
is induced. The exponents P/v, y/u and (1 — ot)/v are obtained at the new fixed point and crossover effects 
are discussed. Surface properties are also studied. 
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1 Introduction 

The study of the influence of impurities on phase transi- 
tions is a quite active field of research, motivated by the 
importance of disorder in real experiments [1,2,3,4]. For a 
disordered system to reach equilibrium, the time evolution 
should be large compared to relaxation processes which 
are themselves governed by the dynamics of impurity re- 
distributions. In practical experiments, such a situation 
can never occur in condensed matter systems and one has 
to deal with quenched disorder [5,6]. 

According to the Harris criterion [7], quenched bond 
randomness is a relevant perturbation at a second-order 
critical point when the specific heat exponent a of the pure 
system is positive. The analogous situation when the pure 
system exhibits a first-order phase transition was studied 
later. Imry and Wortis, generalizing the Harris criterion, 
argued that quenched disorder should soften the transi- 
tion and could even induce a continuous phase transi- 
tion [8]. Phenomenological renormalization-group studies 
inspired from the Imry-Ma argument for random fields [9] , 
suggest that in two dimensions, an infinitesimal amount 
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of randomly distributed quenched impurities changes the 
transition into a second-order one [10,11,12,13,14], while 
in larger space dimensions a finite threshold is necessary 
to produce the same effect. The first exhaustive large- 
scale Monte Carlo study of the effect of disorder at a 
temperature-driven first-order phase transition was per- 
formed by Chen, Ferrenberg, and Landau [15] who studied 
the two-dimensional (2D) eight-state random-bond Potts 
model. This model is known to exhibit, in the pure ver- 
sion, a first-order transition when the number of states 
q is larger than 4 [16], and, the larger the value of q, 
the sharper the transition. This property makes the Potts 
model a good candidate for testing the effect of quenched 
bond disorder. Chen, Ferrenberg, and Landau first showed 
that the transition is softened to a second-order phase 
transition in the presence of bond randomness, and ob- 
tained critical exponents very close to those of the pure 
2D Ising model at the new critical point [17]. Since then, 
different results obtained independently emerged. While 
they confirm the second-order character of the phase tran- 
sition, they conclude to a new universality class [18, 19, 20] . 

The essential properties of random systems are gov- 
erned by disorder fluctuations usually described by nor- 
mally distributed random variables. All physical quanti- 
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ties depend on the configuration of disorder, and the study 
of the influence of randomness requires an average over 
disorder realizations. Among the systems where the pres- 
ence of fluctuations is also of primary importance, aperi- 
odic systems have been of considerable interest since the 
discovery of quasicrystals [21]. Quasipcriodic or aperiodic 
distributions of couplings strengths appear as an alterna- 
tive to quenched bond randomness, albeit built in a deter- 
ministic way, making any configurational average useless. 
Their critical properties have been intensively studied, es- 
pecially in the Ising model (for a review, see, e.g., ref. [22]). 
The characteristic length scale in a critical system is given 
by the correlation length and, as in the Harris criterion for 
random systems, the coupling fluctuations on this scale de- 
termines the critical behaviour. An aperiodic perturbation 
can thus be relevant, marginal or irrelevant, depending on 
the sign of a crossover exponent involving both the correla- 
tion length exponent v of the unperturbed system and the 
wandering exponent u which governs the size-dependence 
of the fluctuations of the aperiodic couplings [23] . In the 
light of this criterion, the results obtained in early papers, 
mainly concentrated on the Fibonacci and the Thue-Morse 
sequences (see e.g. Refs. [24,25,26,27,28]) found a consis- 
tent explanation, since, resulting from the bounded char- 
acter of fluctuations, a critical behaviour which belongs to 
the pure model universality class was found in two dimen- 
sions. 

In the last years, much progress have been made in the 
understanding of the properties of marginal and relevant 
aperiodically perturbed systems. Exact results for the 2D 
layered Ising model and the quantum Ising chain have 
been obtained with irrelevant, marginal and relevant ape- 
riodic perturbations [29,30,31,32]. The critical behaviour 
is in agreement with Luck's criterion, leading to essen- 
tial singularities or first-order surface transition when the 
perturbation is relevant and power laws with continuously 
varying exponents in the marginal situation with logarith- 
mically diverging fluctuations. A strongly anisotropic be- 
haviour has been recognized in this latter situation [33, 
34,35]. 

The effect of quasiperiodic and aperiodic distributions 
of exchange couplings at first-order phase transitions has 
only recently been investigated. It was shown that the 
transition remains first-order for the eight-state Potts model 
on a quasiperiodic tiling [36], while a finite-size scaling 
study using Monte Carlo simulations has shown strong 
evidences in favor of a second-order phase transition for 
the "Paper- Folding" aperiodic perturbation [37]. In the 
present paper, we report an extensive Monte Carlo study 
of the influence of aperiodic modulations of the coupling 
strengths on the nature of the phase transition in the 2D 
eight-state Potts model. We arc interested in both bulk 
and surface properties, and several aperiodic sequences 
are considered. 

The paper is organized as follows: In Section 2, after a 
summary of the essential properties of aperiodic sequences 
and a presentation of the layered structure of the system, 
the critical point of the models is exactly located through 
duality. A qualitative description of the phase transition 



is given in Section 3 from a numerical study of the tem- 
perature dependence of some physical quantities. Eventu- 
ally Section 4 contains the results of a Finite-Size Scaling 
(FSS) analysis. 

2 Layered aperiodic structure and details of 
the Monte Carlo simulations 

The Thue-Morse sequence is an example of aperiodic suc- 
cession of digits fk = or f leading to bounded fluc- 
tuations. It may be defined as a two digits substitution 
sequence which follows from the inflation rule 

0^5(0) = 01, 1^5(1) = 10, (1) 

leading, by iterated application of the rule on the initial 
word 0, to successive words of increasing lengths: {fk} = 
01101001. ..It is well known that most of the 
properties of such a sequence can be characterized by a 
substitution matrix whose elements My are given by the 

number nf^ of occurrences of digits i in the substitution 
S(j) [38]. The largest eigenvalue of the substitution matrix 
is related to the length of the sequence after n iterations, 
L n ~ A™, while the second eigenvalue A 2 governs the be- 
haviour of the _cumulated deviation from the asymptotic 
density p x = f k : 

L 

^(/ fe -Poc)H^r~(^r, (2) 
k=i 

where the wandering exponent is defined by: 



The spin system considered in the following is a layered 
two-dimensional 8-state Potts model. The Hamiltonian of 
the system with aperiodic interactions can be written 

-PH = Y, Ki i 5 w ( 4 ) 

(iJ) 

where the spins ai, located at sites i, can take the values 
a = 1,2, ... ,q, the sum goes over nearest-neighbour pairs, 
and the coupling strengths are allowed to take two differ- 
ent values K = K and K\ = Kr. They are distributed 
according to a layered structure i.e. the distribution is 
translation invariant in one lattice direction, and follows 
the aperiodic modulation {fk} in the other direction: In 
layer k, both horizontal and vertical couplings take the 
same value Kr^ k . This layered structure is reminiscent in 
the shape of the correlated clusters obtained by Monte 
Carlo simulations (Fig. 1). 

Particular choices of coupling distribution make it pos- 
sible to determine exactly the critical point by duality ar- 
guments [37] . Consider a system of N layers with a distri- 
bution {fk}, made from a succession of vertical-horizontal 
(V-H) bonds when read from left to right (Fig. 2), and let 
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Fig. 2. Layered structure of the system and its transformation 
under duality (above). The critical line in the coupling space 
is also shown (below). 

The fluctuations of the coupling strengths per bond at 
the correlation length scale £ induce a thermal perturba- 
tion {St) ~ t-K'"- 1 ), which has to be compared to the 
deviation from the critical point t. The resulting pertur- 
bation has a crossover exponent <P = 1 + v(u> — 1) and 
is relevant when <P > 0. This criterion, first obtained by 
Luck [23], founds its justification at a critical point, i.e. 
when the correlation length of the pure system diverges 
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Table 1. Substitution rules for the aperiodic sequences con- 
sidered in the text. 

Sequence substitutions wandering 

exponent 

Thue-Morse (TM): -> 01, wtm = -oo 

I -> 10. 

Paper-Folding (PF): 00 -> 1000, lo pf = 
01 -f 1001, 
10 -f 1100, 

II -> 1101. 

Three- Folding (TF): -> 010, utp = 

1 -> 011. 



Table 2. Details of the typical parameters used in the 
Finite-Size Scaling Monte Carlo simulations" and energy auto- 
correlation time in the case r = 5. Sizes L — 2 n or 3 n between 
the limits indicated in the table have been used. 

Sequence size (6) # MCS/spin (c) t e 



min. max. 

PS 8 to 256 2.10 b to 3.10 b Tz 651.1 

TM 8 to 2048 1.10 s to 4.10 6 6.7 141.4 

PF 8 to 2048 3.10 5 to 4.10 5 6.7 8.8 

TF 3 to 2187 2.10 5 to 2.5 10 5 4.5 18.2 



as the transition point is approached. The question of its 
application at first-order phase transitions is not yet clear. 
The purpose of this paper is to report some results in this 
situation. The eight-state Potts model has the advantage 
of undergoing a strong first-order transition in the pure 
case, and has been already intensively studied by several 
techniques with random bonds [17,18,19]. 

In the following, we consider three different aperiodic 
sequences and a periodic system (PS) with the regular 
succession of couplings Ki, Kq, K\, Kq, . . ., in which the 
transition is surely first-order. This system constitutes a 
reference for the first-order type behaviour and presents 
the advantage of having the same value for the critical cou- 
pling (at fixed r) than the aperiodic sequences considered. 
The substitution rules for the aperiodic sequences studied 
in the paper are given in Table 1. Details on their proper- 
ties can be found in ref. [29, 33] . We first performed prelim- 
inary runs (to be presented in the next section) for differ- 
ent values of the temperature, and then intensive Monte 
Carlo simulations at criticality on two-dimensional square 
lattices of sizes L x 4L (PF, TM, PS) or L x 3L (TF), us- 
ing the Swendsen-Wang cluster algorithm [41] . This tech- 
nique is known to be very efficient to study second-order 
phase transitions, since it is less affected by the critical 
slowing down than conventional Metropolis algorithm. At 
first-order phase transition points, other methods can be 
used to improve the effectiveness, but we favoured the use 
of the same algorithm to study the different regimes, in- 
creasing the number of Monte Carlo iterations when nec- 
essary in order to obtain reliable results. The multi-spin 
coding technique has also been used to speed up the sim- 
ulations [42]. The geometry L x pL allows a sufficiently 
large number of rows in order to explore the aperiodic 
structure at long enough length scales and the value of 
p has been chosen with respect to the symmetries of the 
sequences. The boundary conditions are periodic in the 
short direction (L) and free or fixed in the long direction 
(pL)} 

With these boundary conditions, translational invari- 
ance holds in the vertical direction, and a local order pa- 



1 The details of the simulations at the critical point (sizes, 
boundary conditions, autocorrelation time, # of MC itera- 
tions) are given in Table 2. 



a 5000 iterations (in MCS/spin) have been discarded. 
6 The values indicated correspond to the largest size pL (hori- 
zontal direction) with p = 4 (PS, TM, PF) and p = 3 (TF). 
The same numbers of MC iterations have been used for the 
two types of boundary conditions (free and fixed in the hori- 
zontal direction). 
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Fig. 3. Temperature dependence of the order parameter profile 
in the case of the Thue-Morse sequence. 

rameter is defined by the majority orientation of the spins 
at column j [45] : 



m(j) 



<lPn 



1 



Ma 



(m(j)). 



(7) 



Here, p max (j) = max (T [p (T (j)], where p a (j) is the density 
of spins in the state a at column j and (...) denotes the 
thermal average over the Monte Carlo iterations. The sys- 
tems under consideration arc highly inhomogeneous, as it 
can be seen in Fig. 3, so, in order to reduce fluctuations, 
we studied average quantities, e.g. 



- 1 



q-1 



M= (to), 



(8) 



where p max has the same meaning as above, over the whole 
system and is not restricted to a given row. The suscepti- 
bility is obtained as usually via the fluctuations of magne- 
tization x = KpL 2 ((m 2 ) — (m) 2 ), and we also computed 
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the energy density 



(9) 



where the prefactor ensures a normalization to 1. Local 
properties at the surface have also been calculated, e.g. 
Mr = (m(l)>. 



3 Off-critical point behaviour 

In this section, we give a qualitative description of the or- 
der of the phase transition. For this purpose, we performed 
preliminary MC simulations over a range of values of K 
for a system of size L x 2L (L from 8 to 256, TM and PF 
sequences), and determined the temperature dependence 
of the physical quantities. The behaviour of the average 
magnetization, susceptibility, and Binder cumulant of the 
magnetization for the TM and PF sequences at r = 5 are 
shown in Fig. 4 and Fig. 5, respectively. 
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Fig. 4. Temperature dependence of the average magnetiza- 
tion, susceptibility, and magnetization cumulant for the TM 
sequence (r = 5). The different symbols correspond to simula- 
tions of systems of sizes 8 x 16 (o) to 256 x 512 (x). 




Fig. 5. Temperature dependence of the average magnetization, 
susceptibility, and magnetization cumulant for the PF sequence 
(r = 5). The different symbols correspond to simulations of 
systems of sizes 8 x 16 (o) to 256 x 512 (x). 



As well known in Monte Carlo simulations, it is dif- 
ficult to observe, in the numerical data, a jump of the 
order parameter at the transition point of a first-order 
phase transition, and similarly the <5-like behaviour of the 
susceptibility cannot easily be distinguished from a pure 
power-law, so we report here also the results for the mag- 
netization cumulant [46] . As a consequence of the highly 
inhomogeneous systems under consideration, these cumu- 
lants exhibit a quite complicated structure, but their be- 
haviour already gives an idea of the nature of the tran- 
sition for the two sequences. One can indeed observe in 
Fig. 4 (TM) that a narrow well appears in the vicinity 
of the transition point and becomes deeper as the system 
size increases. This should be the signature of a first-order 
phase transition, while in Fig. 5 (PF) there is no analogous 
significant trend. The direct comparison between the two 
sequences also shows that the variation of the magnetiza- 
tion and of the susceptibility close to the critical coupling 
is sharper for TM than PF. 

A criterion to analyse the order of phase transitions, 
from the observation of the way the non-analytic behaviour 
develops as the critical point is approached, has recently 
been proposed [47]. Singular quantities, like the suscepti- 
bility, are response functions with diverging non-analytic 



6 



C. Chatelain, P.E. Berche, and B. Berche: Aperiodic Eight-State Potts Models 



2000 




-0.1 0.0 0.1-0.2 0.0 0.2 



K-K(L) K-K(L) 

c c 

Fig. 6. Rescaled susceptibility with respect to the size- 
dependent critical coupling K C (L) for TM and PF sequences 
at r = 5. Three sizes (64 x 128, 128 x 256 and 256 x 512) have 
been used. In the case of TM, one observes a behaviour which 
approaches a 5-function, while a power-law is obtained for PF. 



towards the values and 2. This is characteristic of a first- 
order phase transition, since the scaling dimensions asso- 
ciated to the temperature and magnetic field, y t and ijh, 
respectively, take a special value equal to the dimension 
d of the system [48]. In the case of the PF sequence, the 
behaviour is drastically different, and these effective ex- 
ponents evolve towards non trivial values around 0.5 and 
1. 



4 Finite-Size Scaling 

4.1 Dynamical exponent 

The conjectures of the previous section have to be con- 
firmed by a FSS analysis. As well known in Monte Carlo 
simulations, the energy auto-correlation time te is a good 
test to know about the order of the transition. In Table 2, 
we have given the characteristic values of te for different 
sequences and sizes at the critical point. The number of 
MC iterations is always of order 10 4 te to ensure reliable 
results. The autocorrelation time is shown in Fig. 7 in a 
semi-logarithmic scale. 



behaviour in the thermodynamic limit. In the scaling re- 
gion, there exists a certain interval where such functions 
are decreasing at first-order transitions, leading to cross- 
ings of the rescaled curves which evolve towards a <5-like 
behaviour, while they increase with the size of the system 
in the neighbourhood of K c at second-order transitions 
(power-law behaviour). This is illustrated in Fig. 6 where 
the case of Thue-Morse sequence belongs to the first situ- 
ation, while Paper-Folding corresponds to the second one. 
For this latter perturbation, typical Monte Carlo simula- 
tions are shown in Fig. 1. In the disordered phase, the 
layered structure of the correlated clusters becomes ap- 
pearent as the system approaches the transition point, 
where these clusters begin to grow in the perpendicular 
direction, leading to a second order critical point in the 
thermodynamic limit. 

Similar qualitative observations were reported in Ref. [37] 
where temperature-dependent effective exponents, for av- 
erage magnetization and susceptibility, were computed by 
comparing the data at two different sizes L and L' = L/2: 
In the case of the magnetization for example, assuming a 
scaling form M L (t) = L~ t3 / V M(Lt v ), where t =\K - K c \ 
and M(x) is a scaling function, the quantity 



x L (t) 



In M L /M L 



InL/L 1 ' 

expanded in powers of Lt v close to K c leads to 



X L (t) * - 



f3 | W M'{Lt v ) 
2 In 2 M(Lt v ) 



+ 0{LH lv ) 



(10) 



(11) 



which defines an effective exponent which evolves towards 
— f as the critical point is approached and in the ther- 
modynamic limit. In the case of the TM sequence, the 
successive estimates of § = d — yu and ^ = 2yh — d evolve 
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Fig. 7. Energy autocorrelation time r at K c (r = 5). For TM, 
the dashed line is a fit to an exponential behaviour, while it is a 
power- law fit for PF (z ~ -0.04) and TF (z ~ 0.12) sequences. 
The data corresponding to the periodic system have not been 
fitted. 



The numerical data for PF and TF sequences can be 
fitted by a power-law te ~ L z , with a very small dy- 
namical exponent presumably linked to a logarithmic be- 
haviour, while in the case of TM, te is exponentially di- 
verging te ~ 

L d/2 e 2*L^ where 

a is an order-disorder 
interface tension. These results support strong evidences 
that in the case of Paper-Folding and Three-Folding se- 
quences, the fluctuations are strong enough to soften the 
transition to a second-order regime. In the case of Thue- 
Morse, even if the autocorrelation time, compared to the 
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periodic case, is lowered by the fluctuations, the transition 
remains first-order. We note that the data corresponding 
to the PS were not fitted since one has to go to strongly 
first order transitions (q = 15) to observe the predicted 
exponential behaviour in the pure system [49]. 

4.2 Bulk properties 

We can now enter upon a more refined characterization of 
the phase transition in the two regimes. Since the critical 
point is exactly known, FSS techniques are well indicated 
to get accurate results. We made different series of simu- 
lations with free- free (/|/), fixed- fixed (F\F), and mixed 
(f\F) or (F\f) boundary conditions (BC) in the horizontal 
direction. From all simulations we can extract the singular 
behaviour of the magnetization and of the susceptibility, 
since there is no regular contribution and, at the critical 
point we can write 

M{K C ,L)=A M L-< 3 ' V , (12) 

X (K C ,L) = A X L^, (13) 

where Ax are non-universal critical amplitudes. On the 
other hand, the four series of simulations are necessary in 
order to extract the singularity associated to the energy 
density which contains a regular part including both a 
bulk and a surface x i" 1 contribution [33,50]. 

This latter part must be split in two terms, since the two 
surfaces are different. We thus have 2 

E (Mr) (K c ,L) = EW(K c ,L) 

+ (^ ( - 1) (ir c ,L)+ J B(- 1 )(^ c ,i))xL- 1 
+ A Wr) L^- 1 ^ + ... (14) 

where (l\r) specifies the BC's (free or Fixed) for the left 
and right surfaces. We note that our simulations being per- 
formed in a cylinder geometry, the Euler number vanishes, 
and thus there is no InL term in the free energy [51,52]. 
The asymptotic values of the bulk energy density with dif- 
ferent BC's are the same in the thermodynamic limit, but 
the amplitudes of the finite-size corrections being differ- 
ent, the regular contributions cancel in the combination: 

AE(K C , L) = E (F \ F) + E U \ f) - Ey\ F) - E (F \ f) 

= [Af\f) + A/I/) - A (f\F) - A {F \ f) ]L {a - 1)/v 

(15) 

leaving a pure power-law. 

The scaling dimensions can be deduced from log-log 
plots of the different quantities vs. the system size. In the 
case of the Thue-Morse sequence, a crossover appears, as 
the size increases, towards a behaviour which resembles 
the periodic one, characterized by a vanishing exponent for 
the magnetization. This effect is visible in Fig. 8 where the 

2 We note the the singular surface terms being less divergent, 
they only add a correction to scaling. 



evolution of size-dependent effective exponents is shown. 
Similarly, after the crossover regime, the susceptibility is 
described by an exponent j/v close to the value d = 2. 
The behaviour of the energy density difference exhibits 
log-periodic oscillations, characteristic of systems with dis- 
crete scale invariance. It makes more difficult a quantita- 
tive analysis, 3 but a tendency to a decreasing slope for the 
largest size is nevertheless observed. 

Certainly, a careful analysis is needed to avoid the 
crossover effects from small sizes to the true fixed point 
behaviour in the infinite lattice size limit, so we define 
the following procedure: From the log-log curves between 
pL mm and pL max , one determines an effective exponent 
x(L m i n ) for each quantity; then the smaller size is canceled 
from the data and the whole procedure is repeated until 
only the three or four largest sizes remain. The effective 
exponent is then plotted against £~ in . This prescription 
makes appearent the crossover effects and enables us to 
identify unambiguously the asymptotic regime 4 . Different 
values of the aperiodic perturbation (r = 0.2, 5, and 10) 
are shown in Fig. 8 in the case of TM sequence. At small 
sizes, the system is still strongly under the influence of the 
fluctuations induced by the aperiodic distribution of cou- 
plings, while as the size increases, the effective exponents 
converge towards trivial values which are characteristic of 
a first-order regime. The ratio (3 /is = is indeed charac- 
teristic of a discontinuity of the order parameter, while 
7/1/ — 2 is consistent with this discontinuity and with the 
scaling law 2j3/v + 7 /v = d. This behaviour of effective 
exponents is the signature that the aperiodic fluctuations 
are eventually irrelevant. This is corroborated by the fact 
that the crossover takes place at larger sizes when the 
perturbation amplitude becomes stronger (r = 10). 

On the contrary, the two other aperiodic sequences 
exhibit power-law behaviours with non-trivial exponents. 
Since a second-order phase transition occurs for these se- 
quences, the question of the stability of the new fixed point 
has to be considered. A numerical study at different ra- 
tios of interactions 5 r — 0.2, 2, 3, 5, and 10 shows that the 
corresponding exponents remain stable for strong enough 
perturbations (r not too close to the pure system value 1). 

The numerical results of power-law fits in the linear 
part of log-log plots are given in Table 3 for PS, TM, PF 
and TF sequences. Since the values deduced from these 
power-law fits seem to remain stable with respect to the 
perturbation amplitude, we can use the effective expo- 
nents (Fig. 9) to determine a more accurate value of the 
critical exponents deduced from the extrapolation at infi- 
nite size. The results are given in Table 4 for all sequences 
at r — 5, for which value we have the more exhaustive nu- 
merical results. Log-periodic oscillations again appear in 

3 The discrete rescaling factors for the different sequences 
take the values: 4 (TM), 2 (PF), and 3 (TF). 

4 One can nevertheless mention that the smallest strip size is 
not of great value, since it is smaller than the correlation length 
of the pure model at the transition point in the disordered 
phase [53]. 

5 The simulations are more complete in the case r = 5, since 
we have one size less for the other values. 
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Fig. 8. Effective size-dependent exponents associated to the 
susceptibility and magnetization in the case of TM sequence 
for different values r = K\ /Ko = 0.2, 5, and 10, of the coupling 
ratio (Note that r = 0.2 corresponds to a perturbation of the 
same strength than r = 5). The error bars correspond to the 
standard deviations of the corresponding power-law fits. 



the behaviour of the energy density combination in the TF 
sequence, but the tendency is coherent with the behaviour 
of PF sequence. 

From these values of the critical exponents, we can 
deduce the scaling dimensions associated to the temper- 
ature and the magnetic field, yt — d — « 1.00 and 

Dh = d — & = d+ ^/ v i=y 1.50 at the new fixed point for 
PF and TF. It could be surprising to obtain, within the 
precision of our results, the same fixed point for these ape- 
riodic perturbations, but we can mention here that both 
of them have the same wandering exponent u — 0. 



4.3 Surface properties 

The surface properties can also be investigated, and here 
we ask if the aperiodic perturbations are also liable to 
modify the surface critical behaviour. We determined nu- 
merically the value of the order parameter at both surfaces 
j = 1 and j = pL for the four sequences considered. It 
gives the corresponding exponents, called Pi/v and /3 p l / v, 
respectively. The log-log plots are shown in Fig. 10 and the 
exponents, deduced from the slopes of log- log plots in the 
linear regime, are given in Table 5. 

We can observe that the PS and TM systems exhibit 
analogous singularities at the surfaces. It is known in the 
pure case that, even with a first-order phase transition in 
the bulk of the system, a second-order phase transition is 
obtained at the surface. The critical exponents keep con- 
stant values around 0.6 independently of the interaction 
ratio r. In the case of PF and TF sequences, a second-order 



Table 3. Bulk critical exponents obtained by the slope of 
finite-size scaling results for the fours sequences. The numbers 
in parentheses give the estimated uncertainty in the last digit. 



PS 



TM 





p/v 


■y/v 


1 

P/v 


j/v 


r = 0.2 
r = 5. 
r = 10. 


0.046(6) 
0.072(4) 
0.045(2) 


2.01(1) 
1.97(1) 
2.02(1) 


0.10(3) 
0.08(2) 
cross." 


1 81(7) 
1.90(5) 
cross. 




PF 


TF 






p/v 


■y/v 


P/v 


■y/v 


r = 2. 
r = 3. 
r = 5. 
r = 10. 
r = 0.2 


0.464(7) 

0.475(8) 

0.480(3) 

0.49(1) 

0.46(1) 


1.13(2) 

1.07(2) 

1.017(9) 

1.01(1) 

1.04(2) 


_b 

0.43(2) 
0.44(2) 
0.43(1) 


1.15(2) 
1.09(3) 
1.15(2) 






(1- 


a)/v 






PS 


TM 


PF 


TF 


r = 5. 


0.05(4) 


osc. c 


1.015(4) 


osc. 



a "cross." means that the crossover is still too strong to allow 
any linear regime in the log-log plots. 

b The symbol - means that the corresponding runs have not 
been performed. 

c "osc." means that the log-periodic oscillating behaviour does 
not allow any precise estimation of the exponent. 

Table 4. Bulk critical exponents obtained by extrapolation at 
infinite size of finite-size scaling results for the sequences PS, 
TM, PF and TF with a ratio r = 5. 





PS 


TM 


PF 


TF 


P/v 


0.012 


0.020 


0.499 


0.508 


l/v 


1.986 


1.993 


0.995 


1.009 


(l-a)/v 


0.05 


osc." 


1.001 


osc. 



a "osc." means that the log-periodic oscillating behaviour does 
not allow any precise estimation of the exponent. 



regime which depends on the coupling ratio is obtained. 
The transition is strenghtened when the coupling are suf- 
ficiently enhanced, and one can even observe a first-order 
transition (in the case of PF) at the surface. 5 This re- 
sult can be understood by the behaviour of the average 
coupling at a length scale n in the vicinity of the left sur- 
face for example, K n , compared to the asymptotic average 
coupling Kaoi 



Rn — 



Kn 



Pn (r- 1) + 1 
Pco (r-l) + l ; 



1 " 

On = ~ y^fk- 

k=l 



(16) 



5 This first-order transition is only possible here because the 
system being two-dimensional, the surface alone cannot order 
above the critical temperature where the bulk is not ordered. 
In higher dimensions, this regime should lead to a surface tran- 
sition. 
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y/v 



PS 

TM 

TF 

- - ^ 

PF 



1 



-p/v 



-.,ti- 



(a-l)/v 



0.0 



0.1 



0.2 



7/L 



Fig. 9. Effective size-dependent exponents associated to the 
susceptibility, magnetization, and energy density difference 
(dotted line: PS, dot-dashed line: TM, long-dashed line: TF, 
and solid line: PF) for r = 5. Arrows on the right of the figure 
demarcate the curves corresponding to j/v and —j3/v. 



Table 5. Surface critical exponents" obtained by the slope of 
finite-size scaling results for the four sequences. 







PS 




TM 




fc/v 








r — 5. 
r = 10. 
r = 0.2 


0.60(1) 
0.60(1) 
0.56(1) 


0.59(1) 
0.60(1) 
0.61(1) 


0.55(2) 
cross. c 
0.60(2) 


osc." 

cross. 

osc. 






PF 




TF 








d 





3. 
5. 
10. 



0.0023(6) 
0.0000(0) 
0.0000(0) 



0.58(1) 
0.58(1) 
0.58(1) 



0.530(4) 0.17(1) 
0.532(2) 0.095(3) 



r = 0.2 0.54(1) 0.000(0) 0.19(1) 0.61(2) 



a The numbers in parentheses give the estimated uncertainty 
in the last digit. 

b "osc." means that the log-periodic oscillating behaviour does 
not allow any precise estimation of the exponent. 
c "cross." means that the crossover is still too strong to allow 
any linear regime in the log-log plots. 

The symbol - means that the corresponding runs have not 
been performed. 



10 J 



10 



10' 



10" 



10" 



10"' 



10' 



: --» s . m v. 










M , 


a 


PF 












TM 
























| PS 






M 4L 







10' 



10° 



a few rows (from 2 to 10) in the vicinity of the surface 
reproduces the same exponents. 




pL 



Fig. 10. Local surface magnetization at both ends of the sys- 



Fig. 11. Average coupling at a length scale n from the left 
surface for the sequences PF, TM and TF (r = 2). The inter- 
actions are enhanced, on average, close to the surface for PF, 



tern for the four sequences considered in the text. The curves . ., ., . ', . . . , . ,, 

while they are reduced tor 11 , and remain constant m the case 



have been shifted for clarity. 



of TM. 



This ratio is always greater than 1 when r > 1 for PF 
and produces a significative enhancement of the interac- 
tions close to the boundary, leading to a decrease of the 
exponent of the surface magnetization. The contrary hap- 
pens with TF as shown in Fig. 11. We have also checked 
that these exponents are characteristic not only of the lo- 
cal boundary behaviour, but also of an average surface 
property, since the average of the order parameter over 



4.4 Crossover effect or marginal variation of the 
exponents for PF and TF sequences? 

The question one poses in this section is whether the small 
observed variations of the critical exponents in the second- 
order phase transition regime results from a crossover ef- 
fect or from a marginal behaviour. It is interesting here to 
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Table 6. Bulk and surface critical exponents 2 obtained by the 
slope of finite-size scaling results for PF sequences in the case 
of the Ising model q = 2. 



r 


p/v 


7/1/ 




Pal/v 


0.2 


0.337(2) 


1.37(1) 


0.497(1) 


0.0203(4) 


5 


0.360(2) 


1.351(9) 


0.0176(7) 


0.517(3) 


10 


0.407(1) 


1.233(6) 


0.0044(2) 


0.525(4) 



a The numbers in parentheses give the estimated uncertainty 
in the last digit. 

make a comparison with what occurs in the Ising model 
case q = 2. For this model, the PF and TF sequences are 
known to lead to a marginal behaviour with continuoulsy 
varying critical exponents. The surface properties of this 
system have been intensively investigated [33,35], but the 
coupling distribution used here being different from the 
one used in previous works, we cannot directly compare 
the values of the exponents. We performed a few simula- 
tions for q = 2, and the results are given in Table 6. It 
appears clearly that the values of the different critical ex- 
ponents exhibit, as expected, a continuous variation with 
the amplitude of the interactions. The aim of these sup- 
plementary simulations is to show unambiguously that a 
marginal variation of the exponents is strong enough to 
be distinguished numerically from the crossover effect at 
small perturbations. Our results suggest that the small 
variation observed in the Potts model case is probably 
due to this latter situation. 



5 Conclusion 

In this paper, we have shown that the deterministic fluctu- 
ations generated by an aperiodic modulation of nearest- 
neighbour couplings in the eight-state Potts model pro- 
duce a softening of the transition, and are even liable to 
induce a second-order phase transition. It happens when 
the fluctuations around the average coupling are strong 
enough, and it is the case for the Paper-Folding and Three- 
Folding sequences, although they are characterized by a 
vanishing wandering exponent. In the case of the Thue- 
Morsc sequence, the wandering exponent being —00, the 
transition remains of first-order, and the scaling dimen- 
sions keep their pure values. These results are consistent 
with Luck's criterion [23], provided that we replace the 
correlation- length exponent v by its trivial value l/yt = 
l/d at the first-order fixed point. The crossover exponent 
associated to the aperiodic distribution <P = 1 + (u> — l)/d 
is then positive (relevant perturbation) for PF and TF se- 
quences while it is negative (irrelevant) in the case of TM 
sequence. 

The analysis of the bulk properties shows that the new 
fixed point exponents (PF and TF) are stable, i.e. do not 
depend, up to small crossover effects, on the value of the 
perturbation amplitude. This is clearly different from the 
marginal behaviour encountered for similar sequences in 
the Ising model case. We can furthermore notice that the 



new universality class seems to be robust, i.e. the same 
for both sequences, a result which is not a priori obvi- 
ous. One can nevertheless mention that our results are 
coherent with the stability of the new fixed point, which 
requires a non-positive value for the crossover exponent 
4>' in the second-order regime. With uo = 0, takes the 
value <P' — 1 — v at the new fixed point. Using hypcrscaling 
relation a = 2 — dv (which should hold, unless anisotropic 
behaviour is found), the value a ~ 0, obtained numerically 
for PF, leads to v ~ 1 and thus <P' ~ 0. It seems reason- 
able to propose a renormalization group sketch, illustrated 
by the evolution of the effective exponents, where the two 
possible sets of exponents should correspond to two dif- 
ferent fixed points, the stability of which depend on the 
strength of fluctuations (i.e. the value of the wandering ex- 
ponent). The two types of situations are shown in Fig 12. 
The periodic system, as it can be seen in Fig 9, is not 
influenced by the existence of the second fixed point. 




r=0 



Fig. 12. Evolution of the size-dependent effecive exponents for 
the two types of aperiodic sequences. 

The surface magnetization has also been analyzed and 
analogous conclusions can be given. It is interesting to 
notice that a first-order surface transition (with a second- 
order regime in the bulk) can be induced in the case of 
PF, i.e. the exact contrary of the pure model behaviour. 

There are still some open questions concerning aperi- 
odic perturbations in these systems. A possible anisotropic 
scaling behaviour could be obtained in the second-order 
induced regime, since it occurs in the Ising model already. 
One should also investigate aperiodic sequences with di- 
verging fluctuations characterized by a positive wandering 
exponent, like Rudin-Shapiro for example. A second-order 
transition should also be obtained, but its universality 
class could be different. 
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